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Abstract. Clusters of galaxies contain a fair sample of the universal baryonic mass fraction. A combined analysis 
of the intracluster medium (ICM) within their hydrostatic regions, as derived from both Sunyaev-Zeldovich effect 
(SZE) measurements and X-ray images, makes it possible to constrain the cosmological parameters. We consider 
both gas fraction estimates and angular diameter distance measurements. Adopting median statistics, we find, at 
the 2-a level, the pressureless matter density, ^mo, to be between 0.30 and 0.40 and the Hubble constant, Hq, 
between 44 and 66 Km s~^ Mpc~^. 
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1. Introduction 

Rich clusters of galaxies are the largest known virialized 
structures in the Universe. Their importance in observa- 
tional cosmology is known from early times. In the be- 
ginning of the last century, direct estimates of their to- 
tal masses first stated the need for unseen dark matter 
(Zwickv 11933(1 . Now, they continue to provide important 
information to characterize the Universe (see, for example, 
Wang et allOTI Sereno lMt . 

A method, independent of the nature of the dark mat- 
ter and based on minimal assumptions, to constrain the 
geometry of the Universe and its matter and energy con- 
tent is based on gas mass observations in clusters of galax- 
ies. 

As shown with X-ray and optical imaging, the main 
contribution to the baryonic budget in clusters of galaxies 
is provided by the hot gaseous intracluster medium (ICM) . 
The gas mass is about an order of magnitude larger than 
the mass in stars in cluster galaxies and provides a rea- 
sonable estimate of the cluster's baryonic mass (Fukugita 
et alUSi Lin et aL OnUHjl . 

The ICM is X-ray emitting via thermal bremsstrahlung 
and produces a spectral distortion of the cosmic microwave 
background radiation (CMBR), known as Sunyaev- 
Zeldovich effect (SZE) (Sunyaev & Zeldovich [HFtHI 
Birkinshaw I1999|l . The SZE is proportional to the pres- 
sure integrated along the line of sight, i.e. to the first 
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power of the gas density. X-ray emission depends on the 
second power of the density. The ICM mass fraction may 
be estimated from either of these, fundamentally different, 
observables, based on independent techniques. 

Since there is no efficient way to change the baryon 
fraction averaged within a Mpc scale, the mass compo- 
sition in clusters, when estimated out to a standard hy- 
drostatic radius, is expected to reflect the universal mass 
composition (White et al. lDMI Sasaki lTMBI Evrard lTMTjl . 
Precise measurements of the visible baryonic mass frac- 
tion in galaxy clusters, /b, along with our knowledge of 
the universal baryonic density parameter, rieoi provides 
a physically based technique for estimating cosmological 
parameters (White et al. I1993|l . Ratio of the ICM mass 
to the total mass of the cluster represents a first estimate 
of /b. So, a useful constraint on the cosmological parame- 
ters is given by the identity f2Mo = ^bo/ Ib- Independent 
estimates of /b can be obtained either from SZE measure- 
ments or X-ray imaging observations. The angular diam- 
eter distance, where the dependence on the cosmological 
parameters appears, to the cluster enters in the two gas 
mass estimates through a characteristic length-scale of the 
cluster along the line of sight. From the different depen- 
dencies of SZE and X-ray emission on the density of ICM, 
different dependencies on the assumed cosmology are im- 
plied. 

Galaxy cluster gas mass fractions from SZE measure- 
ments have been used in Grego et al. H2001|l to constrain 
JImo- X-ray gas mass fraction in relaxed clusters have pro- 
vided upper limit on the cosmological density parameter 
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(Mohr et al. lT^ Ettori & Fabian [TO^ Allen et al.[ 
Erdogdu et al. 12002 Ettori et a\. fIUU% . 

Without referring to the universal baryonic fraction, 
with some assumptions about the geometry of the cluster, 
a joint analysis of SZE measurements with X-ray imag- 
ing observations, since the different density dependencies, 
makes it possible to determine the distance to the clus- 
ter. Such a distance, independent of the extragalactic dis- 
tance ladder, is then used to measure the Hubble constant 
(Birkinshaw [T^M Mason et al. IMTTI Jones et al. IW)^ 
Reese et al. 1^017^ . 

All of the methods discussed above are based on ICM 
mass observations. Estimates of gas mass depend on the 
cosmological parameters through the angular diameter 
distance to the cluster. Equating the gas fraction in clus- 
ter to the universal baryonic fraction allows to investigate 
both the Hubble constant and the cosmological density 
parameters. Taking advantage of the different density de- 
pendencies, SZE and X-ray observations provide indepen- 
dent constraints in the space of cosmological parameters, 
leading one to solve for two unknowns at the same time. 

Distance measurements are obtained by equating the 
central densities as derived from SZE and X-ray meth- 
ods. Instead of equating the gas fraction to an universal 
value, now the two gas mass estimates are forced to coin- 
cide each other. This constraint is independent and nearly 
orthogonal to the previous ones and allows to solve for an 
additional unknown. 

Usually, when deriving cosmological constraints from 
gas mass fractions, in order to estimate 57mo, a- prior on the 
Hubble constant is assumed. On the other hand, using the 
cosmic distance scale from interferometric measurements 
of the SZE to determine Hq requires, with the current 
data sets, to fix the background cosmology. We show how 
a joint analysis of all the information on the ICM from 
both SZE and X-ray measurements enables to determine, 
at the same time, JImo and Hq. 

The quality of the data samples used in the analysis 
determines the area of the overlapping region and the pre- 
cision of the estimate. We perform a combined analysis of 
data samples, available in literature, of SZE measurements 
and X-ray observations. 

Since we have to combine different data sets, obtained 
with independent methods, we adopt median statistics. As 
shown in Gott et al. 1)20011 see also Avelino et al. '2002' and 
Chen & R,atra Ennn| . median statistics provide a powerful 
alternative to likelihood methods with fewer assump- 
tions about the data. Statistical errors are not required 
to be known and Gaussianly distributed. Since errors, as 
reported in X-ray and SZE literature, are usually asym- 
metric, performing an analysis without using the errors 
themselves turns out to be a very conservative approach. 
Furthermore, median statistics is also less vulnerable to 
the presence of bad data and outliers. Since we are propos- 
ing a method to constrain cosmology with ICM mass esti- 
mates which extend and combine methods previously es- 
tablished, in this early stage we think that median statis- 
tics can be the best choice. 



In Section |21 we derive the dependence of the ICM 
gas mass, derived from X-ray observations or SZE mea- 
surements, on the assumed cosmological parameters. In 
Section 3, the data samples and our selection criteria on 
the clusters are presented. In Section 4, the median statis- 
tics is shortly presented and the results are listed. Section 
5 is devoted to the discussion of various systematic uncer- 
tainties. Conclusions are in Section 6. 



2. The gas mass estimate 

The ICM mass is calculated by integrating the ICM den- 
sity profile, picM, over an assumed shape. 



Micm{V) = I picM{r)d'^r (X mp / nc(r")d r, 
V Jv 



(1) 



where rie is the electron number density profile and mp is 
the proton mass. For a spherically symmetric system, like 
the widely used /3- model (Cavafiere & Fusco lWHl I1978|l . 
the mass within a fixed metric radius r is 



MicM oc Ucor , 



(2) 



where rieo is the central number density. The radius r can 
be expressed as the product of an observable angular ra- 
dius by the angular diameter distance to the cluster, cIa- 
The central number density can be measured by fitting 
the X-ray surface brightness profile, Sx oc J nlA{Tc)dl, 
where the integration is along the line of sight and A(Te) is 
the X-ray emissivity at the electronic temperature Tg . The 
dependence on distance is made explicit by expressing the 
integration variable in a dimensionless form, dl = dxdQ. 
It comes out 




(3) 



so, for the gas mass, when measured with X-ray observa- 
tions, it is 



(4) 



Spatially resolved measurements of the SZE can also 
determine rioo. Since the brightness temperature decre- 
ment of the CMBR towards a cluster is expressed as 
- / ncTod/, then 



T'CMBR 



AT. 



oc 



sz 



1 



SO, it is 



sz 

ICM 



(5) 



(6) 



The temperatures generally taken in analyses are the X- 
ray emission weighted temperatures, usually measured 
over a few core radii. If the temperature has a spatial 
structure, the temperature inferred from such an average 
procedure may depend on how much of the cluster is con- 
sidered. This can lead to a substantial change in the SZE 
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inferred, and thus to a systematic error in the determina- 
tion of the cosmological parameters (Majumdar & Nath 

Typically, an estimate of the total mass of a cluster 
of galaxies is obtained under the assumption that the 
gas, supported solely by thermal pressure, is in hydro- 
static equilibrium in the cluster's gravitational potential. 
Assuming spherical symmetry and isothermal gas, the to- 
tal mass of a cluster within radius r is 



TOT 



(<r) 



kBTxr dlog nc(r) 
Gfirrip dlogr 



oc dA, 



(7) 



where fee is the Boltzmann's constant, fimp is the mean 
molecular weight of the gas and Tx is the spatially aver- 
aged X-ray emission temperature of the gas determined 
from a broad-beam spectroscopic instrument. Combining 
Eqs. (HEl) with Eq. Q, the two estimates of the ICM 
mass fraction turn out 



ICM 



MicM ,13 

M 

MtOT 



(8) 



et al. 11999(1 : H{z) is the redshift dependent Hubble con- 
stant, 

H{z) = HoVQuoil + z)-^ + f^Ao + f^Ko(l + ■ (11) 

The density contrast dc is attained at the radius rg^ 
(Evrard et al. 119961 Evrard [T997|l . According to an analy- 
sis of gas velocity moments (Evrard et al. I1996|l . Sc = 500 
is a conservative estimate of the boundary between the in- 
ner, nearly hydrostatic and virialized central region of the 
cluster and the surrounding, recently accreting outer enve- 
lope. Clusters of different temperatures have similar struc- 
tures once scaled to rsoo, within which the easily visible 
region of the cluster is also probed with the typical back- 
ground sensitivity. Hydrodynamical simulations (Evrard 
et al. I1996|l also show that the assumption of an isother- 
mal gas in hydrostatic equilibrium is valid within the virial 
radius. 

The virial equilibrium expectations at a fixed density 
contrast, that is 



where, for X-ray observations, (3x = 3/2 and, for SZE 
data, /3sz = 1- 

By eliminating rico from Eqs. (jSJEl), one can solve for 
the angular diameter distance, yielding 



T ' 



dA OC — ^. 



(9) 



The derived values of /icm depend on the cosmolog- 
ical parameters through the angular diameter distances 
(Sasaki I1996|l . In a Friedmann-Lcmaitre- Robertson- 
Walker universe filled in with prcssureless matter and a 
cosmological constant, the angular diameter distance to a 
source at redshift z is 



dAiz) 



xSi 



1 



Ho |f^Kor/2(l + z) 



O |l/2 




where JIao is the reduced cosmological constant, Uko — 
1 — J7mo ~ ^^AOi and Sinn(a;) is sinh(x), x, sin(a;) for fl^Q 
greater than, equal to and less than zero, respectively; Hq 
is the today Hubble constant. For the expression of the 
distance in inhomogeneous universes we refer to Sereno et 
al. I|2()()lll2nn^ . In the next, we will consider only the flat 
case, 17ko = 0, strongly supported by the bulk of evidences 
(de Bernardis et al. 2000 ; Harun-or-Rashid & Roos .2001,1 . 

To compare the ICM mass fraction of different clus- 
ters, we have to study the same portion of the virial re- 
gion in each cluster. In such a region, we expect all clusters 
to have the same gas fraction and the gas to be isother- 
mal. Regions of different clusters are characterized by the 
same properties if they encompass the same mean inte- 
rior density contrast, 6c, with respect to the critical den- 
sity at their own redshit, pc{z) = ^"^^g ; as strongly sug- 
gested by numerical simulation (Evrard et al. [T596i Frenk 



GM ScPcT^ 



oc 



(x6cH^iz)rl, 



(12) 



can be calibrated by means of numerical simulation 
(Evrard et al. .l996j) . 



''50o(T'x) = rio(500,f2Mo) 



lOKcV H^{z) 



1/2 



1 



(13) 



where h is Hq in units of 100 km s^^Mpc^^. The normal- 
ization rio in Eq. H13|) is the average radial scale of 10 KeV 
clusters at density contrast Sc = 500. Since the scaling law 
reflects the virial equilibrium within 6c ~ 0(10^), rio(500) 
is nearly insensitive to the peculiar cosmology. We put 
rio(500,rJMo) = 1.24 ±0.08 Mpc (Evrard et al. [TM5|l . 

The gas fraction at ^500, once known at a flxed phys- 
ical radius rx, can be evaluated using a mild, power law 
extrapolation of the data quoted at rx (Evrard Il997fl . 



/iCM[r5oo(7x)] = /icM(r-x) 



(14) 



with 77 = 0.17 as derived from numerical simulations 
(Evrard et al. ll996|l . The small value of 77 also reduces the 
error on the gas fraction which propagates from an uncer- 
tain determination of the normalization factor in Eq. ((13(1 . 
Using a different normalization based on observations of 
the relatively relaxed cluster A1795, Mohr et al. 1(1999(1 
found rio(500) w 1.19. Since t] = 0.17, the relative vari- 
ation on /icm (500) between the two normalization is less 
then 1%. 

Present X-ray observations are mostly sensitive to the 
inner cluster regions, but, future SZE data would probe 
more exterior region, allowing an analysis on larger scales. 
There are some evidences that the gas fraction increases 
by about 15% from r^oo to 7-200 (Ettori & Fabian [T999(l . 
Since the currently derived gas fraction proflle increases 
regularly from the inner part to the outer part, up to 
the virial radius and beyond (Sadat & Blanchard 12001(1 . 
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the 77 value at Sc "C 500 accordingly changes but, when- 
ever the ratio of rg^ and rx is within few percent of each 
other, the effect of 77 on the final result is really negligible 

(CooraylMi- 

The variations in the gas fraction within the virial 
regions are then quite small. At a given temperature, 
the scatter in /icm is less than 20%, including both in- 
trinsic variations and measurements errors (Arnaud & 
EvrardPgg; Vikhhnin ITOMjl . 

The estimated gas mass fraction out to ^500 depends 
on cosmology through the angular diameter distance, see 
Eq. (|5J), and the time dependent Hubble constant, see 
Eq. (O. From Eqs. (0113111, we get 



/iCM[r5oo] « di-^H{zr^ 



(15) 



/icm[?'50o] is a decreasing function of both ^mo and Hq. 
So, an increment in 77 determines an overestimate of these 
cosmological parameters. 

Under the assumption that the gas fraction is time 
independent and constant for all clusters, it is possible 
to constrain the cosmological parameters (Sasaki 199Q 
Pen [TM7I Cooray I1998al Danos & Pen Rines et 

al. I1999|l . The evolution of the gas mass fraction is still 
under discussion (Schindlcr 1999; Ettori & Fabian .199 9 
Roussel et al. I2000|l . An observed variation of /icm with 
the redshift would be explained by a wrong assumption 
for the angular diameter distance and the cosmological pa- 
rameters and/or a true time evolution. In the near future, 
the study of the power spectrum of the secondary CMBR 
anisotropics due to the thermal SZE by clusters of galaxies 
should give a discriminatory signature of any possible evo- 
lution of the ICM mass fraction fMaiumdar I200l|l : how- 
ever, at the moment, numerical simulations do not suggest 
any evolution. 

3. Data sample 

The observational sample we consider is based on the X- 
ray data set in Mohr et al. (|1999. hereafter MME) and 
Ettori & Fabian ifT^ EE), on the SZE data from Grego 
et al. 1)20011 GCR), and on the joint analysis of interfero- 
metric SZE observations with X-ray imaging observations 
in Reese et al. 1^717^ RCJ). All the published ICM mass 
fractions are given out to an angular radius within which 
the signal-to-noise ratio is good enough and the problems 
deriving from extrapolation procedure are really negligible 
(Sadat & Blanchard'5DM ). 

In order to work with a homogeneous sample, we im- 
pose a cut on the temperature, requiring Tx > 5 KeV. 
Cluster evolution is not an entirely self similar process 
driven exclusively by gravitational instability (David et 
al. lT^ Ponman et al. lTMSI MME). While massive clus- 
ters are less affected by processes like galaxy feedback, low 
mass clusters may have lost g result of preheating 

and post-collapse energy input, so enhancing ICM deple- 
tion within the virial region. The net result is that the ICM 
mass fraction shows a mild increasing trend with the tem- 



perature. Both observations (MME) and numerical sim- 
ulations of the effect of preheating (Bialek et al. I2001|l 
show that /b is depressed below the cosmic mean baryon 
fraction in clusters with Tx ~ 3 KeV. However, the trend 
between /b and Tx is not clear (Roussel et al. 12000(1 and 
an analysis in Arnaud & Evrard H1999|l of selected clus- 
ters with weak cooling flow does not give clear statistical 
significance. Sanderson et al (2003 ) found departures from 
a self-similar behaviour in the scaling properties of a large 
sample of virialized haloes. Both the relation between the 
gas density slope parameter and the temperature and the 
gas fraction data reveal a flattening of the gas density pro- 
files in small sized haloes, consistent with energy injection 
into the ICM by non-gravitational means. A clear trend 
for cooler system to have a smaller gas fraction emerged, 
although, above 4 KeV, the significance of this correlation 
is weakened. 

To be more conservative with respect to these short- 
comings, we consider only clusters with Tx > 5 KeV, 
which are less affected by feedback from galaxy formation 
and where /icm appears to be constant without opposite 
claims. This criterion is passed by 28 clusters from the 
sample in MME, 35 from the sample in EE, 18 from both 
the sample in GRC and RCJ 

4. Data analysis 

Median statistics provide a powerful tool to experimen- 
tal data analysis. Few assumptions about the data and 
their errors are required. Usual statistics assume that i) 
experimental data are statistically independent; ii) there 
are no systematic effects; in) experimental errors are 
Gaussianly distributed; iv) the standard deviation of these 
errors is known. On the other hand, median statistics as- 
sume only hypotheses i) and ii). 

To compute the likelihood of a particular set of cosmo- 
logical parameters, we count how many data points are 
above or below each cosmological model prediction and 
compute the binomial likelihoods. Given a binomial dis- 
tribution, if we perform N measurements, the probability 
of obtaining k of them above the median is given by 



P(fc) = 



-N 



k\{N -k)\' 



(16) 



We perform such a test on the baryonic gas fraction mea- 
sured with X-ray data from MME and EF, on the gas mass 
fraction from SZE observations from GCR and, finally, on 
the cosmological distances measurements in RCJ 

4.1. Baryonic mass fractions 

We count the number of baryon mass fractions that are 
too heavy with respect to (i.e. greater than) the ratio 
f^Bo/^^Mo- ICM is not the only contribution to the cluster 
baryon budget. Other contributions come from the lumi- 
nous mass in galaxies, intergalactic stars and a hypothet- 
ical baryonic dark matter. While the latter two terms can 
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be neglected (Ettori et al. l2U(Jl|l . the typical stellar contri- 
bution to the baryonic mass is between 5 and 20% (White 
et al. [TM^ Fukugita et al. 1998", Roussel et al. POOO). 
The ICM is the most extended mass component in clus- 
ters, while the galaxies are the most centrally concentrated 
one (David et al. 11995(1 . There are some indications of a 
trend of increasing gas mass fraction and decreasing mass 
in optically luminous matter with increasing mass of the 
clusters (David et al. 199.^ David IT?)95|I . but this claim is 
stih debated (Roussel et al. 12117111)1 . Lin et al. piH^I found 
that the total baryon fraction is an increasing function 
of the cluster mass. Estimating the stellar mass using K- 
band luminosity, they also showed how the ICM to stellar 
mass ratio nearly doubles from low- to high-mass clusters. 
However, our choice to consider only massive clusters with 
Tx > 5 KeV gives some cautions in considering the to- 
tal mass in galaxies as a fixed fraction of the cluster gas. 
We have adopted the estimate in Fukugita et al. 1)1998(1 . 

/B = (i+[o.i8lg:i°] hy^)f,cM. 

A determination for Qbo is required. One of the most 
precise determination of the physical baryonic density is 
derived in O'Meara et al. 1)2001 1 By combining measure- 
ments of the primeval abundance of deuterium, as in- 
ferred from high-redshift Lya systems, and theoretical 
predictions of the big-bang nucleosynthesis, they found 
riBoh'^ = 0.0205±0.0018. Measurements of the angular 
power spectrum of the CMBR also provide an estimate 
of the baryon abundance. For the recent WMAP data 
(Spergel et al. 2003), ^Boh^ = 0.024±0.001, depend- 
ing primarily on the ratio of the first to second peak 
heights. When combined with other finer CMBR exper- 
iments and 2dFGRS measurements, WMAP data gives 
^Boh'^ = 0.022±0.001 (Spergel et al. 12(1(1^ . in remarkable 
agreement with the result from Lya forest data. CMBR 
predictions depend on a multi-dimensional fitting proce- 
dure involving all the cosmological parameters together, 
included the parameters we want to determine, i.e. Hq and 
fluo, so we prefer to adopt the estimate from O'Meara et 

al. imnn). 

The ICM mass fraction in Eq. ((15(1 is a decreasing func- 
tion of both Hq and J^moj so, increasing the observed value 
of the baryon abundance, i.e. J^boi determines an under- 
estimate of both Hq and flmo- 

We use the ICM mass fractions to a fixed metric ra- 
dius as reported in Table 4 in MME, Table 2 in EE and 
Table 4 in GCR. Then, we extrapolate to the virial radius 
assuming the reference cosmological model {flmo = 1 a-nd 
h = 0.5 in MME and in EF and f^Mo = 0.3 and ft. = 1 in 
GCR); the prediction of a generic cosmological model is 
obtained according to Eq. ((15(1 . We join the X-ray derived 
mass fraction lists in MME and EF in a single data sample 
with 63 data points. The 1-cr and region 2-a in the flMQ-h 
plane are located by cosmological pairs whose predictions 
are between 28 and 35 times (68.65%^ confidence region) 

^ This and the following probabilities are determined by 
adding the single binomial likelihoods of k overestimates in 
a data sample of A*' entries, according to Eq. (Ilfcil . The 1-cr and 
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0.2 0.4 0.6 0.8 1 
Q.MO 

Fig. 1. Concordance region (white) at the 1-a level. 

Different gray-scaled patches show constraints from dif- 
ferent data samples. Regions labelled with RCJ, GRC 
and MME-I-EF denote constraints from cosmological dis- 
tances, baryonic mass fractions from SZE measurements 
and baryonic mass fractions from X-ray data, respectively. 

or between 24 and 39 times (95.70%), respectively, above 

The SZE derived mass fraction data sample provides 
18 entries. Now, the l-cr and 2-a regions in the flMo-h 
plane are located by points with a number of overesti- 
mates between 7 and 11 (76.2%) and 5 and 13 (96.9%), 
respectively. Since the different dependence on the cosmo- 
logical parameters, the results from X-ray and SZE mea- 
surements locate independent regions in the QMQ-h plane, 
see Fig. (QJ. In any case the points in the SZE sample are 
not enough to significantly constrain cosmological param- 
eters from estimates of gas fraction alone. 

4.2. Cosmological distances 

As seen in Section|51 when combining X-ray and SZE data, 
a third constrain can be obtained without referring to any 
value of Hbo ■ Reese et al. (|2002l) determined the distances 
to 18 galaxy clusters with redshift ranging from z '--^ 0.14 
to 0.78. Now, for each cosmological model, we count the 
number of clusters that are too distant, i.e the number of 
cluster whose distance calculated according to Eq. l)TUI) for 
a set of cosmological parameters is greater than the mea- 
sured value. The confidence region in the f2M0-^ plane 
are determined with the same criteria as for the SZE de- 

the 2-a confidence region are chosen to be symmetric around 
the median and such that the total probability inside them is 
just larger than 68.3% and 95.4%, respectively. 
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0.2 0.4 0.6 0.8 1 
Q.MO 



Fig. 2. Concordance region (white) at the 2-cr level. 
Different gray-scaled patches show constraints from dif- 
ferent data samples. Regions labelled with RCJ, GRC 
and MME-t-EF denote constraints from cosmological dis- 
tances, baryonic mass fractions from SZE measurements 
and baryonic mass fractions from X-ray data, respectively. 



0.31 - nuo ~ 0.34 and 0.55 - ft. - 0.62. At the 2-cr level. 
Fig. ©, we get 0.30 ~ Quo ~ 0.40 and 0.44 ~ /i ~ 0.66. 

5. Systematic effects 

Several systematic uncertainties are involved when deriv- 
ing cosmological parameters from the cluster gas mass ob- 
servations. 

5.1. I CM depletion 

Shocks during cluster formation can drive an ICM deple- 
tion within rsoo . Numerical simulations (Frenk et al. 
and observations (Sanderson et al. 2003) show that the 
gas distribution is more extended than the dark matter 
one. The extended gas structure implies a weakly rising 
baryon fraction with radius (EvrardJ997: EF), /icm ~ 
see Eq. (|14ll . and this increasing trend is stronger in less 
massive clusters (Schindler ll999|l . The ratio T of the en- 
closed baryon fraction to the universal cosmic value within 
a fixed density contrast, 

\"M0 /b(c'c)/ 

has been calibrated by simulations. A modest overall bary- 
onic diminution, about 10%, has been found fEvrard fl997l 
Frenk et al. HM^ . 



rived mass fraction. Despite the numbers of entries for 
the cosmological distances and for the baryonic mass frac- 
tion from SZE data are the same, the former does better 
since the difference dependence on cosmological parame- 
ters. The test on the cosmological distances, very sensi- 
tive to the Hubble constant, provides a nearly orthogonal 
constraint to the method of the baryonic fraction. At the 
l-cr level, we get 0.50 ^ h ^ 0.66; the 2-ct confidence 
range is 0.38 ~ /i ~ 0.70. Using a statistics, where the 
statistical uncertainties have been obtained combining in 
quadrature and then averaging not Gaussian, asymmet- 
ric errors, Reese et al. I|2002|l found, as 68.3% confidence 
range, 0.51 ~ /i ~ 0.64, nearly as constraining as our re- 
sult. 

4.3. Concordance analysis 

To evaluate the allowed range of f^Mo and h, we have to 
consider together the three independent constraints. We 
apply a concordance analysis, only retaining the pairs of 
cosmological parameters which lie within l-cr or 2-cr of 
each individual constraint (Wang et al. 2000). This proce- 
dure is conservative with respect to possibly not well con- 
trolled systematic errors and puts in evidence the most 
effective constraints in delimiting the allowed range. As 
can be seen from Figs. (0121), constraints from cosmolog- 
ical baryonic fraction are nearly orthogonal to that from 
cosmological distances. At the l-cr level, Fig. (Q), we get 



5.2. Density clumping 

Small-scale density fluctuations on X-ray measurements 
of the ICM mass arising from accretion events and major 
mergers can introduce a bias in the determination of /fcM- 
Density clumping on a scale less than the resolution of 
the images causes an enhancement of the X-ray brightness 
by a factor Cn — (rig) / {ue}'^ with respect to a uniform 
smooth atmosphere. Since the surface brightness profile is 
proportional to the square of the central gas density, see 
Eq. |j3Jl, /fcM is overestimated by C^^. Numerical hydro- 

1/2 

simulations show that Cn = 1.16 ± 0.01 (Mathiesen et 
al. 11999(1 . However, the spread in this value, on a cluster- 
by-cluster basis, can be much larger than this. 

Since the value of (ric) is not changed by density 
clumping, the estimate of /icm based on the SZE is not 
affected by this systematic effect. 

Currently, there is no observational evidence of signif- 
icant clumping in galaxy clusters (RCJ). 

5.3. Other effects 

Several other systematic effects can affect the gas mass 
fraction measurements. A magnetic field in the cluster 
plasma might support a non thermal component in the 
X-ray emission. It can be effective in the core regions of 
the cluster but it is unlikely that it plays a major role 
out to the hydrostatic radius fCoorav I1998a|l . The non- 
isothermality of the ICM has also been taken into account 



M. Sereno: Qmo and Ho from SZE and X-ray observations 



7 



(Ettori et al. IMTTl Puy et al. 12009 • The presence of a 
temperature gradient increases the gas fraction estimated 
at the virial radius r5oo, whereas correcting by the con- 
tamination of a magnetic field reduces /icm (Ettori et 
al.[2001). 

Heating and cooUng process may also act in galaxy 
clusters, but they do not change in a significant way the 
estimate of the mean /icm in a numerous enough clus- 
ter sample fCoorav I1908a|l . Cooling processes can alter 
the ICM profiles. Since the SZE depends essentially on 
the pressure profile, a cooling flow can lead to an un- 
derestimation of the cosmological distance (Majumdar & 
Nath l2000|l . Even after excluding 80% of the coohng- 
flow region from the analysis, a ^ 10% overestimation of 
Hq may be in order. 

Cluster geometry introduces an important uncertainty 
in SZE- and X-ray-derived quantities CCoorav FiggSbl Puy 
et al. I2000|l . Complicated cluster structure and projec- 
tion effects cannot currently be disentangled. Projection 
effects of aspherical cluster modelled with a spherical ge- 
ometry broaden the distribution of measured gas mass 
fraction and should be corrected by taking into account 
the distribution of ellipticities for the cluster sample 
(Coorav |1998b|l . The effects of asphericity contribute sig- 
nificantly to the distance uncertainty for each cluster, but 
the determination of the Hubble constant from a large 
sample of clusters is not believed to be significantly bi- 
ased (RCJ). 

These various effects do not act all in the same direc- 
tion; furthermore, they require a very detailed modelling 
in order to correct for a not very significant amount. Since 
an ensemble of hydro-dynamical cluster simulations has 
shown that the departures from an isothermal, spherical 
ICM do not introduce serious errors (GRC), we have not 
considered them in our analysis. 

Considering the effects discussed in Section H5.ll 15.2(1 , 
ICM mass fraction should be corrected according to 



/] 



/b — ./icm 1 
and 

-•SZ _ fSZ 
Jb — JICM 

6. Conclusions 



/gal 



1 



1 



ICM 



Jgal 

/icm 



T(500)cy2' 



T(500)' 



(18) 



(19) 



The simple argument presented here, based on the dif- 
ferent dependence of the galaxy cluster ICM mass esti- 
mates on the cosmological parameters as derived either 
from SZE measurements or from X-ray observations, has 
made it possible to infer both the value of the pressureless 
matter density parameter and the Hubble constant. The 
only additional information is the value of the universal 
baryonic matter density. To our knowledge, this is the first 
time that gas mass estimates are used to derive flmo and 
Hq together. 

We have followed median statistics. This technique as- 
sumes only that the measurements are independent and 



free of systematic errors. Physical quantities in litera- 
ture, such as X-ray emission temperature ICM mass frac- 
tions and cosmological distances, are often presented with 
asymmetric errors bars. Supposing these errors following a 
Gaussian distribution with known standard deviation, as 
required by usual analysis, is a hard hypothesis. On the 
other hand, median statistics, based on fewer assumptions, 
can provide very interesting and constraining results, as 
shown in several astrophysical cases (Gott et al. 120011 
Avehno et al. IW)^ Chen & Ratra 12003 ■ To minimize 
the role of systematic uncertainties, we have performed a 
concordance analysis of the data (Wang et al. l20U0|l . 

Our estimates agree with the currently favoured model 
of Universe (Wang et al. l2U0U|l . derived from observational 
constraints such as measurements of the anisotropy of the 
CMBR (Spergel et al. 12003(1 , large-scale structure obser- 
vations (Peacock et al. 12001(1 and evidences coming from 
type la supernovae (Riess et al. 119981 Pcrlmutter et al. 



We have found a value of Qmo in full agreement with 
the bulk of evidences (Harun-or-Rashid & Roos 120011 
Chen & Ratra 01)1)3(1 . Median statistics analyses of var- 
ious collections of measurements has been used in Chen 
& Ratra 1(2003(1 to determine an estimate of I^mo- They 
found 0.20 ~ r^Mo ~ 0.35 at two standard deviations, in 
full agreement with our estimate 0.30 ^ fluo ~ 0.40. 

Our determination of the Hubble constant is indepen- 
dent of the local extragalactic distance scale. We find 
0.44 ~ /i - 0.66 at 2-(T, in agreement with the estimate 
of the Hubble Space Telescope Key Project, h = 0.72±8 
(Freedman et al. 12001(1 . Together with SZE-derived dis- 
tances, time delays produced by lensing of quasars by fore- 
ground galaxies also provide a tool to determine Hq inde- 
pendent of the extragalactic distance ladder. This tech- 
nique suffers by very uncertain systematics but gives re- 
sults in agreement with our result (Witt et al. 12000(1 . We 
remark how both these two global methods tend to yield 
smaller estimates of Hq than the determination by CMBR 
measurements: WMAP data give h — 0.72 ± 0.05 (Spergel 
et al. 12003(1 . The value of the Hubble constant, as de- 
rived applying median statistics to a collection of nearly all 
available pre-mid-1999 estimates, is h = 0.67±0.02 (95% 
statistical) ±0.05 (95% systematic) (Gott et al. 1200111 . in 
agreement with our result. 

In the near future, by increasing the SZE data sam- 
ple, the estimates of the cosmological parameters obtained 
with our method should greatly improve. 
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